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^ ■ The host galaxy of the quasar SDSS J114816.64+525150.3 (at redshift 

^ . z=6.42, when the Universe was <1 billion years old) has an infrared luminos- 
^ ! ity of 2.2x10^^ Lq^'^, presumably significantly powered by a massive burst of 

^ I star formation'^'^'^'^. In local examples of extremely luminous galaxies such 
^ ' as Arp 220, the burst of star formation is concentrated in the relatively 
Q ■ small central region of < 100 pc radius^'^. It is unknown on which scales 
0\ . stars are forming in active galaxies in the early Universe, which are likely 
^ ! undergoing their initial burst of star formation. We do know that at some 

I early point structures comparable to the spheroidal bulge of the Milky Way 
^ ' must have formed. Here we report a spatially resolved image of [CII] emis- 
■ sion of the host galaxy of J114816. 64+525150. 3 that demonstrates that its 
star forming gas is distributed over a radius of ~ 750 pc around the centre. 
The surface density of the star formation rate averaged over this region is 
~1000 Moyear^kpc^. This surface density is comparable to the peak in 
Arp 220, though ~2 orders of magnitudes larger in area. This vigorous star 
forming event will likely give rise to a massive spheroidal component in this 
system. 

The forbidden ^P3/2 fine-structure line of ionized Carbon ([CII]) at 158 

microns provides effective cooling in regions where atomic transitions cannot be excited, 
and therefore helps gas clouds to contract and form stars. [CII] emission is thus known 
to be a fundamental diagnostic tool of the starforming interstellar medium^'^''. Given 
the very bright continuum emission of the central accreting black hole of quasars in 
optical and near-infrared wavebands, standard star formation tracers (such as hydrogen 
recombination lines) cannot be used to study star formation in these systems. The [CII] 
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line is however much brighter than the underlying far-infrared (FIR) continuum, thus 
making it a prime choice to characterize star formation in quasar host galaxies. 

We used the IRAM Plateau de Bure interferometer to resolve the [Cll] emission 
from the z=6.42 host galaxy of Jl 14816.64+525150.3 (one of the most distant quasars 
known^^'-^^; hereafter: J1148+5251) with a linear resolution of ~ 1.5 kpc. J1148+5251 
is one of only two sources for which the detection of [CII] emission is reported at high 
redshift to date^'^^. A large reservoir of molecular gas (2 x 10^'^ Mq), the prerequisite 
for star formation, has been characterized in this system through redshiftcd rotational 
transition lines of carbon monoxide (CO)^'^'^^. At a redshift of z=6.42 the age of the 
Universe was just ~ 870 miUion years (or l/16th of its present age) and 1" on the sky 
corresponds to 5.6 kpc^^'^^. 

The distribution of the [CII] emission is shown in the middle panel of Figure 1. 
Gaussian fitting to the spatially resolved [CII] emission gives an intrinsic source size 
of 0.27"±0.05" (1.5±0.3 kpc). The [CII] emission is embedded within the molecular 
gas reservoir traced by CO^^, however the [CII] emission is offset to the north from 
the optical quasar and the CO peak by ~0.1" (~600 pc). Given the good agreement 
between the position of the optical quasar and the simultaneous 158 micron continuum 
observations (Fig. 1, left) we do not attribute this offset to inaccurate astromctry. 
The significance of the [CII] detection is high enough that it shows spatially resolved 
velocity structure (red and blue contours in the right panel of Figure 1). 

The (rest-frame) FIR continuum emission underlying the [CII] fine is detected at 10 
sigma significance in the integrated frequency spectrum (Fig. 2). If the FIR continuum 
was due to the (unresolved) optical quasar, a 10 sigma point source is expected at the 
optical position. However, from Figure 1 (left) we only find ~50% of the fiux to be 
coincident with the optical quasar position. This implies that the sensitivity of our 
observations is not high enough to image the remaining FIR flux that is presumably 
due to the more extended emission from star formation. This would imply that at most 
50% of the FIR emission can be attributed to heating by the central black hole, i.e. the 
FIR emission is significantly powered by star formation (in good agreement with the 
molecular gas^'^, dense gas^^, radio continuum^ and dust properties-*^'^ of this source). 
In the following we thus assume a FIR luminosity due to star formation of ~1.1 xlO^^ 
L0, i.e., a star formation rate of ~1700 Moyr"^ (assuming a standard initial stellar 
mass function^'^^). The low significance of the resolved FIR emission is the reason why 
it cannot be used to constrain the size of the starburst region. 

The compactness of the [CII] emission implies that massive star formation is con- 
centrated in the central region with radius 750 pc of the system, even though molecular 
material is available on larger scales (but our [CII] observations cannot rule out star 
formation at lower surface densities over the entire molecular gas reservoir). Given the 
star formation rate derived above, we find an extreme average star formation rate sur- 
face density of ~1000 Mq yr~^ kpc^^ (~ 7x10^^ Lq pc~^) over this central 750 pc radius 
region. Similarly high starburst surface densities are also found in the centre of local 
ULIRGs such as Arp220 (where each nucleus of size ~ lOOpc has Lfir=3x1O^^L0), al- 
beit on spatial scales that are by two orders of magnitudes smaller^'^. For comparison, 
the Galactic young starforming cluster associated with Orion KL also exhibits such 
high densities in its central region^° (Lfir=1.2x10^ Lq, area: ~ larcmin^, 0.013 pc^, 
resulting in ~ 10^^ Lq kpc~^), however over an area that is 8 orders of magnitudes 
smaller than in J1148-F5251. 
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In the context of other galaxies in the early universe, this kpc-scale 'hyper'- 
starburst has a star formation rate surface density that is one order of magnitude 
higher than what is found in massive starforming z~2.5 sub millimeter galaxies^^. It 
is however consistent with recent theoretical descriptions of (dust opacity) Eddington 
limited star formation of a radiation pressure-supported starburst on kpc scales^^. The 
high star formation rate surface density is also compatible with other theories describ- 
ing 'maximum starbursts'^^: stars can form at a rate limited by SFR=e x Mgas/tdyn, 
where e is the star formation efficiency, Mgas is the gas within radius r and t^yn is the dy- 
namical (or free-fall) time, given by \Jr^/ (2 G M). For r=750pc, M~Mgas ~1O^*^M0 
a star formation efficiency of e ~0.4 is required to explain star formation rate density 
that we observe in the case of J1148-I-5251. Such high efficiencies may be expected 
given the high dense gas fractions found in local ULIRGs^^. In this calculation, we 
assume that the stellar initial mass function in this object is not significantly different 
from what is known locally. Such a high star formation efficiency could be expected if 
J1148-I-5251 were to undergo a major merger, where the gas is tunneled to the central 
1.5 kpc on rapid timescales. We note however that our observations do not provide 
clear evidence for a merging system and that other mechanisms may be responsible 
for fueling the ongoing starburst^^. We also note that the star formation rate surface 
density of ~1000 Mq yr'^ kpc'^ is a value averaged over the central ~kpc, i.e. this 
value could be significantly higher on smaller scales, which in turn may violate the 
theoretical descriptions of 'maximum starbursts'. 

Our observations provide direct evidence for strong, kpc-scale star formation episodes 
at the end of Cosmic Reionization that enable the growth of stellar bulges in quasar 
host galaxies. Such 'hyper starbursts' appear to have an order of magnitude higher 
star formation rate surface densities on kpc scales than previously studied systems at 
high redshift^^ . The observations presented here are currently the best means by which 
to quantify star formation rates and their surface densities in quasars at the earliest 
cosmic epochs. They thus demonstrate that [CIl] observations will play a key role in 
studies of resolved star formation regions in the first Gyr of the Universe using the 
upcoming Atacama Large Millimeter/submillimeter Array (ALMA)^^. 
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Figure 1: [CII] observations of the z=6.42 quasar J1148+5251 obtained with the 
IRAM Plateau de Bure interferometer. Observations were obtained in the most ex- 
tended antenna configuration during three tracks in early 2007 and 2008 (z/obs=256.17 GHz, 
^'rcst= 1900.54 GHz), resulting in a resolution of 0.31" x 0.23" (l.Tkpcx 1.3kpc; the beam- 
size is shown in light blue colour in the middle panel). The resolved CO emission from 
VLA observations^^ is displayed as colour scale in all three panels. The cross indicates 
the absolute position (uncertainty: 0.03") of the (unresolved) optical quasar as derived 
from Hubble Space Telescope observations^^. Left: Contours represent the far-infrared 
continuum emission obtained from the line-free channels of the [CII] observations in- 
tegrated over a 445kms~^ bandwidth (contour levels are -0.9 (grey), 0.9 and 1.8 mJy 
(black); rms noise: 0.45 mJy). There is good agreement between the optical quasar 
and the peak of the continuum emission, as well as the peak of the molecular gas emis- 
sion traced by CO, demonstrating that our astrometry is accurate on scales of < 0.1". 
Middle: Contours show the [CII] emission over a velocity range of -293 to +293kms~^ 
(contours are plotted in steps of 0.72 mJy; rms noise: 0.36 mJy). The (rest-frame) 
beam-averaged peak brightness temperature of the [CII] emission is 9.4±0.9K (from 
the peak flux of 7.0±0.36mJy at a resolution of 0.31" x 0.23"), which is similar to 
the CO brightness temperature (8.3 K)^^. If the intrinsic temperature of the gas were 
similar to that of the dust^ (30-50 K), this would imply that we have not fully resolved 
the CO or the [CII] emission. Right: Contours of blue- and a red-shifted emission 
(averaged over velocities from 75-175 km s~^ on either side) are plotted as blue and red 
contours at 3.2 and 4.8 mJy, respectively (rms noise: 0.63 mJy). The dynamical mass 
of ~ 10^° Mq within the central 1.5 kpc deduced from these observations (assuming 
Vrot ~ 250kms~^) is in agreement with earlier estimates on larger spatial scales^^. 
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Figure 2: Spatially integrated [CII] spectrum of the z=6.42 quasar J1148+5251. The 
[CII] line is detected at high significance (bandwidth covered: 1 GHz, or llOOkms"^) 

and is present on top of a 4.5±0.62mJy continuum (consistent with an earlier estimate^ 
of 5.0±0.6 mJy). Gaussian fitting to the line gives a [Cll] peak fiux of 12.7±1.05 mJy, a 
full width at half maximum (FWHM) velocity of 287±28kms~^ and a central velocity 
of 3±12kms-^ relative to the GO redshift^ of z=6.419 (i/obs=256.17GHz). This leads 
to a [GII] flux of 3.9±0.3 Jykms"^ (consistent with earlier, unresolved observations^ of 
4.1±0.5 Jykms"^), which corresponds to a [Cll] luminosity^''' of Lj(.jj] = 1.90±0. 16x10^° 
Kkms~^pc~^ or L[cii]=4.18±0.35x 10^ (adopting a luminosity distance of Dl=64 
Gpc^^), yielding L[cn]/LFiR=l-9xlO~^. This ratio is by an order of magnitude smaller 
than what is found in local starforming galaxies (a finding consistent with local ultra- 
luminous infrared galaxies, ULIRGs^'^^'^^). The line-free channels of the [GII] obser- 
vations are used to construct a continuum image of J1148-I-5251 at 158 microns (rest 
wavelength) as shown in Fig. 1 (left). 



